Very large dielectric response of thin ferroelectric films with the dead layers 
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We study the dielectric response of ferroelectric (FE) thin films with "dead" dielectric layer at the 
interface with electrodes. The domain structure inevitably forms in the FE film in presence of the 
dead layer. As a result, the effective dielectric constant of the capacitor e^fj increases abruptly when 
the dead layer is thin and, consequently, the pattern of 180-degree domains becomes "soft". We 
compare the exact results for this problem with the description in terms of a popular "capacitor" 
model, which is shown to give qualitatively incorrect results. We relate the present results to fatigue 
observed in thin ferroelectric films. 
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We have shown recently that the dead layer forming 
at the interface between ferroelectric (FE) thin film and 
electrodes has drastic effect on the electric response of a 
capacitor [Q. It has direct bearing on fatigue observed 
in FE capacitors since in many cases the deterioration of 
the switching behavior, like the loss of the coercive force 
and of the squareness of hysteresis loop, were attributed 
to the growth of a "passive layer" at the ferroelectric- 
electrode interface ^-|^. It is of principal importance 
that the presence of a dead layer, no matter how thin 
in comparison with thickness of the ferroelectric layer, 
triggers a formation of the domain structure in FE film 
. We have shown that when the thickness of the dead 
layer d is not very small, the apparent (net) polarization 
Pa of the ferroelectric with 180— degree domain walls fol- 
lows an approximate relation dPa/dE oc e^/d, which is 
in good correspondence with available experimental data 
(see e.g. |6|j^ and references therein). Importantly, the 
response of this structure to an external bias voltage be- 
comes more rigid when d increases, i.e. when the dead 
layer grows, even in the absence of pinning by defects. 

The implication for real systems is that with the 
growth of the passive layer the hysteresis loop very 
quickly deteriorates and looses its squareness, as ob- 
served. The approximate 1/d dependence of the response 
suggests that the effective dielectric constant of the ca- 
pacitor 



may become very large when the layer is thin. Here C 
is the capacitance of the electroded FE film of area A, 
with L the separation between electrodes. Indeed, when 
the dead layer is thin, the domain width a becomes very 
large, it grows exponentially with 1/d^ ||l[. The response 
of this domain structure is very soft, and this should 
translate into very abrupt increase of the dielectric con- 



stant eg// of the capacitor. It is easy to show that in 
the present case (180° domains) the linear response is 
not changed by electromechanical effect, the change in 
Ee// only appearing in quadratic terms in external field, 
but here we are interested in zero-field value of eg/ / only. 
We have assumed a quadratic coupling between the elas- 
tic strains and the polarization (as in perovskites). In 
this case the linear response of 180° domain structure 
without pinning is not affected by intimate contact be- 
tween the ferroelectric and the dead layer (excluding 
mere renormalization of materials constants by homo- 
geneous stresses). Note that this is invalid in the case of 
90° domains Q or a linear coupling between the strain 
and the polarization ^J. 

Here we address the anomalous behavior of the dielec- 
tric constant ee// in detail. We also discuss the important 
issue regarding the relation of the present results to so- 
called "capacitor" model The "capacitor" model 
(incorrectly) presumes that the dielectric response of the 
dielectric layer is not affected by the presence of the dead 
layer, so the system looks like capacitances in series. We 
show that the effective "dielectric constant" of the FE 
layer, e/, as found in the "capacitor" model [^,||J§1, is ac- 
tually negative. In spite of this the system remains stable 
(stiffness of the domain pattern is positive) . 

The reason for this apparently unusual behavior (e/ < 
0) is that the "dielectric constant" of the FE layer e/ 
is the non-local quantity which characterizes the whole 
system. This non-local behavior is due to long range 
Coulomb interaction, which makes the response rigid 
even when the FE itself would have a negative "dielectric 
constant" (cf. Ref. (l^. Fig. 1). The "capacitor" model 
neglects this essentially nonlocal behavior and, therefore, 
is hardly applicable to the problem of dielectric response 
of thin ferroelectric films, and to the problem of fatigue 
for that matter. The present considerations remain valid 
until the period of the domain structure remains much 
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smaller that the lateral size of the film (or the grain size). 
Those will set some cutoff for the effects considered be- 
low. 

We shall find the response of a ferroelectric film under 
a bias voltage U with thickness I separated from the top 
and bottom electrodes by passive layers with thickness 
d/2 (Fig. i, inset) with the use of the thermodynamic 
potential M: 



F = F, + Fes, 



electrodes 



dVE^/{8TT) - ea(/5a 



(2) 
(3) 



Here Fg is the surface energy of the domain walls, Fes 
is the electrostatic energy, E is the electric field, while 
Sa i'Pa) is the charge (potential) on the electrode a. The 
last term in Eq. (H) accounts for the work of the external 
voltage source(s). Note that the results would not change 
qualitatively is there were only one dead layer, since the 
accompanying depolarizing field would have the same ef- 
fect. 

We are interested in a case when the ferroelectric has 
the spontaneous polarization || z (Fig. 1, inset), i.e. 
we are considering the case of 180— degree domain walls. 
The uniaxial ferroelectric film has the dielectric constant 
Ec (fa) in the z-direction (in the xy-plane), and the di- 
electric constant of the passive layer is eg. We select the 
X-axis perpendicular to the domain walls. The poten- 
tial If, which is related to the electric field in usual way, 
E ~ — V(/9, satisfies the following equations of electro- 
statics in the ferroelectric and the passive layer 



9z2 



= 0, 



= 0, 



(4) 



with the boundary conditions Lp = —{+)U/2, z 
+ {-){l + d)/2, and 



dz 



fg-^ = 47rcr(x), (fif = Lpg, at z^ 1/2, (5) 



where the subscript / (g) denotes the ferroelectric (pas- 
sive layer, or a vacuum gap). Here a is the density of 
the bound charge due to spontaneous polarization at the 
ferroelectric-passive layer interface, a = Pns = i^s, de- 
pending on the normal direction of the polarization at 
the interface, alternating from domain to domain. Fig. 1 
(inset). Thus, we assume that the absolute value of the 
spontaneous part of the polarization Ps is constant in all 
domains and only its direction is alternating from domain 
to domain. We have also assumed a usual separation of 
linear and spontaneous polarization, so that the displace- 
ment vector is — eikE^. + AnPgi, where i,k — x, y, z, 
and the dielectric response ejj. in uniaxial, Fig. 1 (inset). 
Within this approximation we obtain 



F,. 



dAuip 
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where a = zLPs is the density of the bound charge given 
by only the spontaneous polarization at the interfaces 
between the FE and the dead layer, and the integration 
goes over these interfaces. The periodic pattern consists 
of c— domains with widths ai and 02, and the period 
T ~ ai + a2- The solution of the equations (^) is then 
readily found by Fourier transformation 

aix) ^^ake"'"', <Pc(a;, z) = ^ (/?fc„(z)e^'=^, (7) 
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[1 — exp (zfcai)] , 7^ 0, 



(Jk=0 =0-0 = Ps ; = PsO, 

ai + 02 
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where k = 27rn/T = 7m/a, n = 0, ±1, . . ., a = (oi -|-a2) /2 
and the index a — f, g marks the quantities for the FE 
layer and the dead layer. Note that by definition Pa = (Tq 
is the net spontaneous polarization of the FE layer. We 
obtain 
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(13) 



eg Anlao + eJJ 
47r e„l + e^d 



(14) 



is the net charge density induced on the electrode at z — 
L/2, with A the area of the film. One can calculate the 
apparent dielectric constant of the whole capacitor from 



(15) 



The term Fh in Eq. ( p^ ) is due to the net polarization of 
the FE film induced by the bias voltage U, with the first 
term in Fh ( pi] ) corresponding to the term A: = 0, singled 
out in the stray field energy Finh (p^). With the use of 
dfe from m) we obtain 
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(6) where Z?„ = -Dfc„, kn = im/a. 
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The total free energy of the domain pattern per unit 
area is 



eJ + €cd 
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where 7 = is the surface energy of the domain waU, 
with A the characteristic microscopic length [|l|. This 
free energy allows one to determine the equilibrium prop- 
erties and response of the domain pattern to external 
field. 

We can determine a linear response of the system to 
bias voltage from the total energy (|7|), which, for small 
bias [/, can be expanded up to terms quadratic in I] and 
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where cto = PsS, Fstrayip) is the usual stray energy, 
Fstrayia)^ Finh{a; ctq = 6 = 0), givcu by the first term 
in Eq^dlq). This expression results from expanding the 
Eqs.([l2[),(|l6|) in powers of U and S (note that ctq = Ps^)- 
Here 5* is the stiffness of the domain structure with re- 
spect to external bias voltage U [we have omitted the 
constant term cx C/^ in Eq. (p^)]. We have 
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(19) 
(20) 
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where Sink is the (negative) contribution of the inhomo- 
geneous (stray) energy, corresponding to the second term 

in 



inh 



16a 



-ir 



n—l 



(22) 



We can now consider the following limiting cases: 

Thick dead layer (a <^ d).— There we can replace both 
coth in (^3|) by unity, with the result 



F{S = 0) 
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with e = y/eaec 
mum energy for the equilibrium (Kittel) domain width 



eg. The domain structure has a mini- 



a = ax 



V28C(3) 
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so that 



^ ^, a usual relation. Then inho- 



mogeneous contribution to stiffness Sinh ' 
Eq.(|2|), and we obtain from Eqs. (1^,(11) 
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When d > a^, one can neglect the second (stray) term 
in this expression for stiffness, thus recovering our ap- 
proximate expression for the net spontaneous polariza- 
tion Pa = ctq of the ferroelectric film from dlq) ||l[ 
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This approximation breaks down for thinner films but it 
is obvious that the response becomes softer for thinner 
dead layers. The breakdown of this approximate behav- 
ior has also been noticed by Kopal et al. for similar model 
pT[ , but they have not analyzed what actually happens 
to the response at very small thicknesses d of the dead 
layer. Note that Smh < 0, Eq. (^, and, consequently, 
the following inequality always holds: Pa/U > eg/And. 

Thin dead layer (o ^ d).— In this case the replacement 
of coth by unity is not allowed. For a ^ I and y/e^^ — eg 
we have Z3„ = (1 + cothkd/2) and 
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where b = nd/a, Lin{z) = Y^^=i^^l^^ HI]- -^Y using 
the asymptotic expansion of Li„(z) we obtain for d <^ a 
our earlier result M 
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which corresponds to very abrupt increase of the domain 
width to values a ^ qk, with ax the Kittel width (|24|), 
when the dead layer is very thin. Fig. 1. This approx- 
imation, however, is not sufficient to make an estimate 
of the stiffness for very thin dead layers. Indeed, for- 
mally the stiffness there becomes small and negative, 
= 4:nd[l/ (eg + ecd/l) - 1/eg] < 0. This fact simply in- 
dicates that for thin dead layers the stiffness diminishes, 
5^0 and has to be calculated more accurately. The 
exact results, illustrating the abrupt softening of the di- 
electric response for small d, are shown in Fig. 1. 

One may wish to interpret the approximate result for 
slope of the net polarization Pa c>c 1 /d, given by Eq. ( p6| ) , 
in terms of the "capacitors in series" model ||,|. In- 
deed, a similar result follows if one were to assume that 
the capacitance of the dead layer, cx l/d, dominates at 
small thicknesses of the layer d. 

An opinion has even been voiced that the effective di- 
electric constant of the FE layer is infinite (c/ = oo?) 
since the domain walls in our model are not pinned . 
However, such an interpretation would be incorrect since 
e/, as found in the "capacitor" model, is not infinite, but 
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FIG. 1. The dielectric constants of the ferroelectric with 
the dead layers of thickness d/2, as defined in the text. W is 
the domain wall thickness. Note an abrupt increase of the ef- 
fective dielectric constant of the capacitor e^ff in comparison 
with ei = tgL/d, and that the "dielectric constant" of the FE 
layer is negative, e/ < 0. The relation te// ~ ei, which follows 
from the "capacitor" model, is clearly violated when the dead 
layer is thin and the domain width a is large (a ^ ax, the 
Kittel width, and ee// becomes ^ ei). Inset shows schematics 
of the electroded ferroelectric film (capacitor) with the dead 
layers. 



finite and actually negative. To establish this, one has to 
find the voltage drops across the dead layer and FE film. 
The homogeneous part of the electric fields inside the 
FE (dead) layer Ef {E^) and the corresponding voltage 
drops are found from the standard equations 



Edd + Efl = U, 

egEd = tcEf +AT:Pa, 



(30) 



where from one can, for instance, easily recover the ex- 
pression (^4|) for the net charge density on electrodes. 
One obtains the "capacitor" model by assuming that 
the interface between the FE film and the dead layer 
is equipotential and can be viewed as the metallic film 
separating two areas. Simple calculation gives 
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where Uf = Efl. Since always Pa/U > Sg/And, Eq. (26), 
we obtain < 0, Fig. 1. In spite of formally negative 
"dielectric constant" of the ferroelectric layer e/, Eq^dSM, 
which is an artifact of the "capacitor" model P,^,|6|,|13[, 
the system remains stable (stiffness of the domain pat- 
tern is positive). Indeed, shifting the domain walls would 



create a net electric field in the capacitor, see Eq. ( |Tl| ) 
and second term in Eq. (18), and its energy is the source 
of finite stiffness of the domain pattern, even when the 
walls are not pinned, as is very well known since 1960 . 

The reason for this apparently unusual behavior (e/ < 
0) is that the "dielectric constant" of the FE layer e/ 
m (|l]) is the non-local quantity which characterizes the 
whole system: it depends on the properties of the dead 
layer. This non-local behavior is due to long range 
Coulomb interaction, which makes the response rigid 
even when the FE itself would have a negative "dielec- 
tric constant" (for analogous situation in FE films with 
depletion charge see Ref. @, Fig. 1). Thus, the "ca- 
pacitor" model actually operates with obscure quantities 
without much physical meaning. 

We illustrate the difference between the exact results 
for eg/ / and that following from the "capacitor" model 
eeff « ei = tgL/d in the assumption e/ = oo in 
Fig. 1. The softness of the domain structure, charac- 
terized by Pa/U, increases very abruptly when the dead 
layer is thin and the width of the domains is a 3> clk 
[Kittel width (p^J. According to (|3^), e^ff increases 
abruptly and becomes ^ ei in this region. Fig. 1, in 
stark deviation from the prediction of the "capacitor" 
model Eeff ~ ei (l3|. For thick "dead layers" the effec- 
tive dielectric constant e^ff becomes comparable to ei, 
as we showed earlier . 

We think that the above clearly demonstrates a danger 
of applying a naive electric circuit analysis to FE systems 
where the addition of one "circuit element" (dead layer) 
radically changes the electric response of the other, FE 
layer, by introducing a domain structure. It is not sur- 
prising, therefore, that such an approach cannot explain a 
fatigue observed in FE films (see e.g. |^,^ and references 
therein). Certainly, there might be various reasons for 
the fatigue in FE capacitors. We simply observe that the 
growth of the dead layer at the interface with electrode 
makes the dielectric response of the film rigid even when 
the domain walls are not pinned. The present mecha- 
nism gives a correct order of magnitude for the tilt of 
the hysteresis loops therefore the growth of passive 
layer might be the main source of fatigue. In the limit- 
ing case of thin dead layers the period of the ferroelectric 
domains in the absence of pinning becomes much larger 
than the standard Kittel width. As a consequence, the 
domain pattern becomes very "soft" in the absence of the 
pinning of the domain walls, and its contribution to the 
dielectric response becomes very large, since the domains 
with polarization parallel to the external field can easily 
grow at the expense of the domains with the opposite 
polarization. 
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